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ABSTRACT 

We discuss the reionization of the intergalactic medium (IGM) by quasar sources at high redshift. We 
compute the integrated UV background from observed QSOs, taking into account the hydrogen opacity 
associated with intervening Lya clouds and Lyman limit systems. We note that the published data 
appear to indicate a significant underdensity of absorption systems in the Lya forest with column densities 
Nhi > 10 15 cm“ 2 . This deficit results in a reduction of the opacity of the universe by a factor of 1.5-3 at 
z — 3 — 5 relative to previous estimates. The QSO contribution to the metagalactic flux at the Lyman edge 
may be as large as 7912 ( 2 ) « 6[(1 + z)/4.5]° 5 x 10 _22 erg cm -2 s _1 Hz -1 sr” 1 for qo = 0, and slightly lower 
for qo = 1/2. For a density of the diffuse component of the IGM of fiz)/i 2 0 < 0.025, QSOs could photoionize 
a smooth IGM sufficiently to satisfy the constraints imposed by the Gunn-Peterson effect. The epoch of 
reionization could be as recent as z 5. As a result, neutral patches of IGM would be detectable in the 
spectra of high redshift quasars. The patches would appear as absorption line systems with typical column 
densities of 10 19 — 10 20 cm -2 , and velocity widths of 100 — 1000 km s" 1 . 


SUMMARY 


In the past few years, the growing number of observed quasars with z > 4 (Schneider, Schmidt, & 
Gunn 1989) has promoted a renewed interest in the topic of the reionization of the universe. The absence 
of an absorption trough on the blue side of the Lya emission line in the spectra of these high redshift QSOs 
requires the diffuse component of the IGM to have been highly ionized by z ~ 5 (Schneider, Schmidt, & 
Gunn 1991a). We examine whether QSOs detected in optical surveys can supply the required metagalactic 
ionizing flux. 

The key quantity here is J 912 , the mean specific intensity of the diffuse radiation field at the observed 
wavelength \ 0 b 3 = 912A, as seen by an observer at redshift z 0 \ )i : 
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where e(v, z) is the proper volume emissivity (ergs cm -3 s -1 Hz -1 ) of QSOs at frequency v = c(l + 
z)/[912A(l + z 0 &«)] and redshift z, Ho = 50/isokms -1 Mpc -1 , and exp[— r e //(912, z 0 b s , z)] is the average 
transmission over all lines of sight of a clumpy medium. We assess the contribution of quasars to the 
ultraviolet extragalactic background at high redshift in light of the two major uncertainties: (1) the existence 
and degree of a controversial decline at the faint end of the QSO luminosity function for z Z, 3 (Boyle 1991; 
Irwin, McMahon, & Hazard 1991; Schmidt, Schneider, &; Gunn 1991; Warren, Hewett, &, Osmer 1991); 
and (2) the continuum opacity of intervening absorption systems (e.g., Sargent, Steidel, fc Boksenberg 
1989; Lanzetta 1991). We use the QSO luminosity function of Boyle (1991), and consider two cases: (1) a 
population of QSOs with a constant comoving density for z > 2 and which turn on at a redshift z on — 6 
(case CC); (2) a QSO population which declines by a factor of 2 in comoving density per unit redshift for 
z > 3 (case ED). We estimate for the density of the diffuse intercloud medium which must be photoionized 
0.02 < Q d < 0-04 for Ho = 50 and 0.005 < &d < 0.01 for Ho = 80, after subtracting observed optically 
thick absorbers from the estimate of the baryonic density Qfe/i 2 0 = 0.05 ± 0.01 from Walker et al. (1991), 
based on nucleosynthesis constraints. A more complete discussion is presented by Meiksin & Madau (1992). 
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We evaluate the expected UV background from QSOs using two models for the amount of attenu- 
ation by intervening absorption systems. The model used by Madau (1992) gives a moderate amount of 
attenuation. After examining the H I column density distribution in the spectra of 4 QSOs observed at high 
resolution, however, we find a significant underdensity of Lya systems with lo gNpi > 15 (Figure 1). A 
deficit by a factor of 3, significant at the 4cr level, occurs in the range 15 < log Nhi < 16- No Lya forest 
lines were detected for log Nhi > 15, with a significance of 30-. We show the best fit line of slope ft = 1.7 
determined over the range 13.75 < log Nhi < 14.8 ( solid line), and the same fit renormalized to systems in 
the range 14.8 < log Nhi < 15 ( dashed line). The deficit, if real, would result in a decrease in the amount 
of attenuation by a factor of 1.5-3 from the moderate case. We note that the deficit may be a consequence 
of a Jeans instability in the Lya clouds. Spherical clouds of diameter D would be gravitationally unstable 
for column densities greater than Nhi, max = 3.0 x 10 16 (Z)/200 kpc)“ 3 cm “ 2 , where 22 is the UV 
background intensity at the Lyman edge in units of 10 -22 ergs cm" 2 s” 1 Hz^ 1 sr” 1 . 
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We show the resulting UV intensity for case CC in Figure 2. The intensity for both the low ( solid 
line), and moderate ( dashed line), attenuation models satisfies the Gunn-Peterson limits on the Lya optical 
depth, top , for O^/i 2 0 < 0.025. The horizontal bars are lower limits to the required intensity based on 
TGp{{zabs) = 2.6) < 0.05 (Steidel & Sargent 1987), and TGp((z a bs) = 3.8) < 0.31 (Jenkins & Ostriker 1991). 
In each pair, the upper (lower) bar corresponds to Ho = 50(80)kms _1 Mpc -1 . We assume Qp^lo = 0.025 
and a temperature of 15,000 K. Case ED succeeds for the low attenuation model, but fails for a moderate 
amount of attenuation unless Ho > 80. 

The Gunn-Peterson limit places a direct constraint on the epoch of reionization, through the 
porosity parameter Q(z) of the source H 11 regions, when radiative recombinations in the IGM can be 
neglected. For case CC, we obtain TGp(z)Q(z)/Clohl 0 = 0.45/i5 O 1 T 4 ^ 0 - 75 [(3'Fa)/as](H-9o)~ 1 ( 1 + 2: ) 4 ” 2?o [ 1- 
(x/x on ) l+qo ]/J 9 i 2 (z) i where x = 1 + z, a and as are the spectral indices of the UV background and the 
sources, respectively, and T4 is the IGM temperature in units of 10 4 K. The quantity J 9 \ 2 is a dimen- 
sionless integral given by dividing the integrand in eq.(l) by e(c/912A, z 0 bs)- To 15% accuracy, Jon{z) ~ 
0.65(H-z)°- 5_9 ° for the low attenuation model, and 3.0 (l + z)“0+fo) for the moderate model. For case ED, we 
obtain, to 15% accuracy for 3 < z < 6, r G p(z)Q(z)/QDhl Q — 0A2h^QT^ 075 [(i-\-a)/a s ](l-\-z) 3 ~ 2qo /J 9 i 2 (z), 
where Joi 2 (z) ~ 0.35(14- z)°- 75_<7 ° for low attenuation, and 1.2(1 4- z )-(0-5+9o) for moderate. The constraint 
tqp < 0.05 at z — 2.6 requires Q(z — 2.6) > 54(25) for qo = 0(0.5) (a = 0.5, as = 0.7, T4 = 1.5) 
for case CC. The porosity evolves according to Q(z) = <3(0)(1 + z)"( 1+9 °)[l — ( x/x on ) 1+q ° ]. The limit on 
porosity constrains the reionization epoch zj, when Q = 1, to z/ > 5. 9(5. 7) for q 0 = 0(0.5). For case ED, 
Q(z = 2.6) > 40(21) for g 0 = 0(0.5). To 15% accuracy, Q(z) = <2(0)(1 + ^-G+tfo) exp[-0.69(£ - 3)], for 
3 < z < 6. We then obtain zj > 6. 2(5. 2) for <70 = 0(0.5). Thus the reionization epoch may have been as recent 
as zj £ 6, and possibly more recent if radiative recombinations are significant. (Note that this conclusion 
generalizes to any class of radiation sources, depending only on their evolution and spectral shape.) For 
QSOs, we find zj = 5. 8(5. 6) for <70 = 0(0.5) for case CC. For case ED, we find z/ = 6. 2(5. 2) for qo = 0(0.5). 

The epoch of reionization will first reveal itself as absorption lines arising from un-ionized patches of 
IGM. For Q(z) £ 1, the patches will have neutral hydrogen column densities of 10 19 — 10 2O cm“ 2 and black 
cores of velocity width 100 — 1000 kms -1 , determined by the Hubble expansion across a patch. For model 
ED, the lines would persist over an extended range of A z ~ 0.5 — 1, until Q falls below ~ 0.5 and the lines 
broaden into the Gunn-Peterson absorption trough shortward of the Lya emission line. 
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